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Abstract 

Neutrino produced in a chain of nuclear reactions in the Sun starting from 
the fusion of two protons, for the hrst time has been detected in a real¬ 
time detector in spectrometric mode. The unique properties of the Borexino 
detector provided an oppurtunity to disentangle pp-neutrino spectrum from 
the background components. A comparison of the total neutrino flux from 
the Sun with Solar luminosity in photons provides a test of the stability of 
the Sun on the 10^ years time scale, and sets a strong limit on the power 
production in the unknown energy sources in the Sun of no more than 4% of 
the total energy production at 90% C.L. 
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1 Introduction 


The solar photon luminosity (a total power radiated in form of photons into 
space) is determined by measuring the total solar irradiance by spacecrafts 
over the wide subrange of the electromagnetic spectrum, from x-rays to radio 
wavelengths; it has been accurately monitored for decades. The luminosity 
Lq = 3.846 X 10^® W is measured for a precision of 0.4% with the largest 
uncertainty of about 0.3% due to disagreements between the measurements 
of different satellite detectors DEI- The energy lost by neutrinos adds Ly = 
0.023 • Lq to this value [3]. The solar luminosity constraint on the solar 
neutrino fluxes can be written as |1]: 


Lq 

47r(la.M.)2 




( 1 ) 


where 1 a.u. is the average earth-sun distance, the coefficient a, is the 
amount of energy provided to the star by nuclear fusion reactions associated 
with each of the important solar neutrino fluxes, <hj. The numerical values 
of the a’s are determined to an accuracy of 10“^ and better. 

The estimated uncertainty in the luminosity of the Sun corresponds to 
less than 3% uncertainty in total solar neutrino flux. 

The Sun is a weakly variable star, its luminosity has short term fluc¬ 
tuations The major fluctuation occurs during the eleven-year solar 

cycle with amplitude of about 0.1%. Long-term solar variability (such as the 
Maunder minimum in the 16th and 17th century) is commonly beleived to 
do not exceed the short term variations. 

Because of the relation ([T]) between the solar photon and neutrino lumi¬ 
nosity, the measurement of the total neutrino luminosity will provide a test 
of the stability of the Sun at the time scale of 40000 years [6] , the time needed 
for the radiation born at the center of the Sun to arrive to its surface. Finding 
a disagreement between Lq and Lj, would have significant long term enviro- 
mental implications, and in the case of an agreement of two measurements it 
would be possible to limit the unknown sources of the solar energy, different 
from the known thermonuclear fusion of light elements in the pp-chain and 
CNO-cycle. 

The main neutrino sources in the Sun are the pp- and ^Be reactions, 
providing roughly 91 and 7% of the total neutrino flux respectively. Borexino 
already measured ^Be neutrino flux with 5% precision [7] , but till recent time 
the pp-neutrino flux was derived in a differential measurement using the data 
of solar detectors. 

Solar pp neutrinos measurement is a critical test of stellar evolution the¬ 
ory, discussion of the physics potential of the pp solar neutrino flux mea- 
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surement can be found in [SI El [TU] (at the moment of the discussion the 
parameters space for MSW solution was not established yet, thus the au¬ 
thors were giving priority to this part of the physical potential). 

A number of projects aiming to perform pp-neutrino detection have been 
put forward in past two decades, but with all the time passed since the 
proposals, none of them started the operation facing the technical problems 
with realization. The principal characteristics of the proposals are presented 
in table [H The radiochemical experiments sensitive to the solar pp- neutrinos 
(SAGE m] and GALLEX [I2]) are not cited in the table, the combined best 
£t of the radiochemical and other solar experiments gives solar pp-neutrino 
flux of (6.0 ± 0.8) X 10^° cm“^s“^ [TT] in a good agreement with expected 
value of 6.0 x (1.000 ± 0.006) x 10^*^ cm“^s“^. 


Project 

(reference) 

Method 

Threshold 

[keV] 

Resolution 

Mass [tj 

Reaction 

PP 

events 

[d-i 1 

LENS 

fT^ 

LS 

301 (i^ ) 

7% 

@ 1 MeV 

20 

(8% nat ) 


0.5 

INDIUM 

mi 

LS 

118{i^ ) 

5-10% 

@1 MeV 

4 

In -h n'e — 
^^^Sn*(613) -h e- 

1.0 

GENIUS 


ll(e- ) 

0.3% 

1 

i' + e~ 

1.8 

fisl 

scatt 

59(i/ ) 

@ 300 keV 

10 

L' -h e~ 

18 

HERON 

fT^nn 

superfluid He 
rotons / phonons+uv 

50(e“ ) 
141(i/ ) 

8.3% 

@364 keV 

10 

i' + e~ 

u e~ 

14 

XMASS 

fTHl 

liquid Xe 
scintill 

50(e- ) 
141(i^ ) 

17.5% 

@ 100 keV 

10 

1 / -h e~ —>• 

u e~ 

14 

CLEAN 

im 

liquid Ne 

20(e- ) 

82(i^ ) 


135 

u + e~ 
i' -h e~ 

7.2 

HELLAZ 

He (5 atm), 

TPC 

100(e“ ) 
217(i^ ) 

6% 

@800 keV 

2000 

1 / -h e~ —>• 

u -f- e~ 

7 

MOON 

drift 

chambers 

168{i^ ) 

12.4% FWHH 
@ 1 MeV 

3.3 

Ve Mo -> 

^°°To + o- 

1.1 

MUNU 

TPC,CF4 

100(e“ ) 

16% FWHH 

0.74 

e~ 

0.5 

[2^ 

direction 

217(i^ ) 

@ 1 MeV 

(200 ) 

u e~ 


NEON 

mi 

He,Ne 

scintill 

20(e- ) 

82(i^ ) 

16% FWHH 
@ 100 keV 

10 

1 / -j- e~ —V 

V -|- e~ 

18 

10 t LS 

LS 

170(e“ ) 
310(i^ ) 

10.5 keV 
@ 200 keV 

10 

1 / -h e~ —>• 

1 / -h e~ 

1.8 

Borexino 

ED 

LS 

165(e“ ) 
305(i^ ) 

5% 

@1 MeV 

75.5 

(fiducial) 

1 / -t- e~ —>• 

V -f- e~ 

13.6 


Table 1: Key characteristics of the solar neutrino projects sensitive to the 
pp-neutrino. The number of expected neutrino is calculated for the fraction 
of the neutrino spectrum above the threshold, but the region of sensitivity 
(limited i.e. by signal to backgrounds ratio of 1) could be stricter. 


A possibility to use ultrapure liquid organic scintillator as a low energy 
solar neutrino detector for a first time was discussed in [2HE5]. The authors 
come to the conclusion that a liquid scintillator detector with an active vol¬ 
ume of 10 tons is a feasible tool to register the solar pp-neutrino if operated 
at the target level of radiopurity for Borexino and good energy resolution 
(5% at 200 keV) is achieved. 
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Borexino Detector 



Figure 1: Schematic view of the Borexino detector. 

2 Borexino detector 

The Borexino detector consists of a dome-like structure (see Fig{T]), 16 meters 
in diameter, Filed with a mass of 2,400 tons of highly purihed water which 
acts as the shield against the external radioactive emissions of the rocks and 
the environment that surround the facility. The water buffer acts also as an 
effective detector of the residual cosmic rays. Within the volume of water 
a steel sphere is mounted which hosts 2,200 looking inward photomultiplier 
tubes providing 34% geometrical coverage. On the outer side of the stainless 
sphere 200 PMTs of the outer muon veto detector are mounted, these PMTs 
detects the Cherenkov light caused by muons passing through. 

The sphere contains one thousand tonnes of pseudocumene. Finally, the 
innermost core of the facility contains roughly 280 tons of the scintillating 
liquid bounded within a 100 /rm thick nylon transparent bag with ~4.2 m 
radius. The water and the pseudocumene buffers, as well as the scintillator 
itself, have a record-low level of radioactive purity. The energy of each event 
is measured using light response of the scintillator, and the position of the in¬ 
teraction is determined using timing information from the PMTs. The latter 
is important for the selection of the innermost cleanest part of the detector 
within 3 meters radius, as only the internal 100 tonnes of scintillator have 
the radioactive background low enough to allow the solar neutrino detection, 
the scintillator layer close to the nylon serves as an active shield against the 
7 originating from the nylon trace radioactive contamination. The threshold 
of the detector is set as low as possible to exclude triggering from the random 
dark count of PMTs. The Borexino has excellent energy resolution for its 
size, this is the result of the high light yield of ~500 p.e./MeV/2000 PMTs. 
The energy resolution is as low as 5% at 1 MeV. 
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3 Data processing 

The low-energy range, namely 165-590 keV, of the Borexino experimental 
spectrum has been recently carefully analyzed with a purpose of the pp- 
neutrino flux measurement ^7\. The data were acquired from January 2012 
to May 2013 and correspond to 408 live days of the data taking. These are 
the data were collected at the beginning of the second phase of Borexino 
which had started after the additional purification of the liquid scintillator 
following the calibration campaign of 2010-2011 [28]. The main backgrounds 
for the solar neutrino studies were significantly reduced in the Phase 2, the 
content of ®®Kr is compatible with zero, and background from ^^‘^Bi reduced 
by a factor 3 to 4 compared to the values observed at the end of the Phase 
1 just before the purification. 

The experimental spectrum is presented in Fig.l The main features of the 
experimental spectrum can be seen in the hgure: the main contribution comes 
from the decays at low energies (below 200 keV), the monoenergetic peak 
corresponds to 5.3 MeV a—particles from ^^‘^Po decay. The statistics in the 
first bins used in the analysis is very high, of the order of 5 x lO^events, 
demanding development of the very precise model for the studies - the al¬ 
lowed systematic precision at low energy part should be comparable to the 
statistical fluctuations of 0.14%. 

3.1 Data analysis 

The Borexino spectrum in the low energy range is composed mainly of the 
events from decays of present in liquid organic oscillator in trace 
quantities, its measured abundance with respect to the is (2.7 ±0.1) x 
10“^® g/g. The (3 -decay of C is an allowed ground-state to ground- 
state (O’*" —)■ ) Gamow-Teller transition with an endpoint energy of Eq = 

156.476 ± 0.004 keV. 

In Borexino the amount of the active PMTs is high (~2000), demanding 
setting of the high acquisition threshold in order to exclude detector trigger¬ 
ing from random coincidence of dark count in PMTs: hardware trigger was 
set at the level of 25 PMTs in coincidence within 30 ns window, providing 
negligible random events count. The acquisition efficiency corresponding to 
25 triggered PMTs is roughly 50% and corresponds to the energy release of 
~50 keV. In present analysis, the same as in the “pp”-analysis, the threshold 
was set at the lowest possible value at ~60 triggered PMTs (~160 keV). In 
independent measurement with laser the trigger inefficiency was found to be 
below 10“® for energies above 120 keV |26j . 
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Figure 2: Borexino energy spectrum between 165 and 590 keV. The pp neu¬ 
trino component is shown in red, the spectrum in dark purple and the 
synthetic pile-up in light purple. The large green peak is ^^‘^Po a-decays. ^Be 
(dark blue), pep and CNO (light blue) solar neutrinos, and ^^°Bi (orange) 
spectra are almost flat in this energy region. 

3.2 Energy resolution 

The most sensitive part of the analysis is the behaviour of the energy resolu¬ 
tion with energy. The variance of the signal is smeared by the dark noise of 
the detector (composed of the dark noise from individual PMTs). In order to 
account for the dark noise the data were sampled every two seconds forcing 
randomly fired triggers. Some additional smearing of the signal is expected 
because of the continuously decreasing number of the PMTs in operation. 
The amount of live PMTs is followed in real time and we know precisely 
its distribution, so in principle this additional smearing can be precisely ac¬ 
counted for. It was found that the following approximation works well in the 
energy region of interest: 
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a% = N{po - pivi) + N'^{vt{N) + Vf), 


where N = Nq < f{t) >t is average number of working PMTs during 
the period of the data taking, f{t) is a function describing the amount of 
working PMTs in time with /(O) = Nq. The last parameter here is Vf{N) =< 
pit) >T — < fit) >y, it is the variance of the f(t) function over the time 
period of the data taking. 

An additional contribution to the variance of the signal was identified, it 
is the intrinsic width of the scintillation response. From the simple consider¬ 
ation this contribution reflects the additional variations due to the fluctua¬ 
tions of the delta- electrons production and the energy scale non-linearity. It 
should scale inversely proportional to the energy loss. Because of the limited 
range of the sensitivity to this contribution, basically restricted to the very 
tail of the spectrum, the precise energy dependence could be neglected 
and we used a constant additional term in the resolution. Taking it all to¬ 
gether, the variance of the energy resolution (in terms of the used energy 
estimator) is: 


a% = N{po - pivi) + N^{vt{N) + Vf) + 

where is contribution of the dark noise (fixed) and atnt is contribution 
from the intrinsic line shape smearing. The probability piis linked to the 
energy estimator with relation n = Npi. 

3.3 Scintillation line shape 

The shape of the scintillation line (i.e. the response of the detector to the 
monoenergetic source uniformly distributed over the detector’s volume) is 
another sensitive component of the analysis. A common approximation with 
a normal distribution is failing to describe the tails of the MC-generated 
monoenergetic response already at the statistics of the order of 10^ events. 
This was notihed already in the first phase of Borexino and the approximation 
of the scintillation line shape with generalized gamma- function ra have been 
used to £t monoenergetic ^^°Po peak in the solar ^Be neutrino analysis m 
[3I1171. The generalized gamma- function (GGC) was developed for the energy 
estimator based on the total collected charge, but it provided a reasonable 
approximation for the energy estimator based on the number of triggered 
PMTs given the moderate statistics corresponding to the total amount of the 
events in ^^°Po peak. The fit quality of the ^^‘^Po peak is rather insensitive 
to the residual deviations in the tails. This is not the case for the precise 
^^G spectrum modeling, as all the events in the fraction of the ^^G spectrum 


above the energy threshold originate from the spectral smearing. An ideal 
detector’s response to the point-like monoenergetic sonrce at the center is a 
perfect binomial distribntion and it would be well approximated by a Poisson 
distribution. When dealing with real response one should adjust the “base 
distribution” width to take into account at least the additional smearing 
of the signal due to the various factors. The problem with binomial “base 
function” (or with its Poisson approximation) is that its width is defined by 
the mean value. In case of Poisson the variance of the signal coincide with 
mean /i. A better approximation of the response function was tested with 
MC model, namely the scaled Poisson (SP) distribution: 

fix) = (2) 

(xsj! 

featuring two parameters, that could be evaluated using expected mean and 
variance of the response: 


s = — and 11 = —. (3) 

n al 

The agreement of the approximation and the detector response function 
was tested with Borexino MC model and it was found that at low energies 
(j2|) better reproduces the scintillation line shape compared to the generalized 
gamma function up to the statistics of 10®, while at energies just above the 
tail both distribution gives comparable approximation. The quality of 
the fit was estimated using criterion, for example with 10^ total statistics 
(these events are uniformly distributed in the detector and then the FV is 
selected) for n = 50 (approximately 140 keV) we found x^/''^-<^-/-=88-0/61 
for the GGC compared to x^/n.(i./.=59.3/61 for the SP distribution. 

As proven by MG tests, the SP distribution as a base function works well 
in the energy region of interest despite of the additional smearing due to the 
factors enlisted in the previous paragraph. This is a result of the “absorption” 
of the relatively narrow non-statistical distributions by the much wider base 
function, as follows from MG such an absorption results in the smearing of 
the total distribution without changing its shape. 

As it was noted above, the £t was performed in n scale. All the theoretical 
spectra involved in the fit were first translated into the n scale and then 
smeared using resolution function (|2]) with /r and scale factor s calculated 
using (|3]). As it is clear from the discussion above, the detector’s response 
has the shape described by ([2]) only in the “natural” n scale. If the measured 
values of n would been converted into the energy, the shape will be deformed 
because of the non-linearity of the energy estimator scale with respect to the 
energy, complicating the construction of the precise energy response. 
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3.4 Standard fit 


The “standard” options of the spectral £t are: nnniber of triggered PMTs 
in a fixed time window of 230 ns (npmts) nsed as energy estimator; 62-220 
npmts fit range; 75.5±1.5 tonnes fidncial volnme (defined by the condition 
R<3.02 m and jZj <1.67 m). 

The rate of the solar nentrinos is constrained either at the valne fonnd 
by Borexino in the different energy range (R(^Be) = 48 ± 2.3 cpd 0) , or 
hxed at the prediction of the SSM in the SMW/LMA oscillation scenario 
(R(pep)=2.80 cpd, R(CNO)=5.36 cpd). All connts here and below are 
qnoted for 100 tonnes of LS. 

The rate was constrained at the valne fonnd in independent mea- 
snrement with the second clnster R(^'^C) = 40 ± 1 Bq (or i?(^^C) = 
(3.456 ± 0.0864) x 10® cpd). The synthetic pile-np rate was constrained 
at the valnes fonnd with the algorithm. The normalization factors for other 
backgronnd components were mainly left free (*^Kr,^^°Bi and ^^°Po) and 
the fixed rate of (i?(^^^Pb) = 0.06 cpd) was calcnlated on the base 

of the amonnt of identihed radon events. The light yield and two energy 
resolntion parameters (vt and (Tint) are left free. The position of the ^^°Po 
is also left free in the analysis, deconpling it from the energy scale. 

3.5 Systematics study 

An evident sonrce of systematics is nncertainty of the FV. The FV mass is 
defined nsing position reconstrnction code, residnal bias in the reconstrncted 
position is possible. The systematic error of the position reconstrnction code 
was defined dnring the calibration campaign, comparing the reconstrncted 
sonrce position with the nominal one [271ES]. At the energies of interest the 
systematic error on the FV mass is 2%. 

The stability and robnstness of the measnred pp nentrino interaction rate 
was verihed by performing hts varying initial conditions, inclnding £t energy 
range, method of pile-np constrnction, and energy estimator. The distribn- 
tion of the central valnes for pp-nentrino interaction rates obtained for all 
these £t conditions was then nsed as an estimate of the maximal systematic 
error (partial correlations between different factors are not exclnded). 

The remaining external backgronnd in the fidncial volnme at energies 
relevant for the pp nentrino stndy is negligible. In the particnlar case of the 
very low-energy part of the spectrnm, the fit was repeated in five smaller 
fidncial volumes (with smaller radial and/or z-cut), which yields very similar 
results, indicating the absence of the influence of the external backgrounds 
at low energies. 
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4 Results and Implications 



Figure 3: Survival probabilities for electron neutrino (Borexino only data 

from [271 El 1321E3]) 

The solar pp neutrino interaction rate measured by Borexino is pp = 
144 ± 13(sfat) ± 10{syst) cpd/100 t, compatible with the expected rate of 
PPtheor = 131 ± 2 cpd/100 t. The corresponding total solar pp-neutrino flux 
is (j)pp{Borex) = (6.6 ± 0.7) x 10^° cm“^s“^, in a good agreement with the 
combined best fit value of the radiochemical and other solar experiments 
(ppp{other) = (6.0 ±0.8) x 10^*^ cm“^s“^ [H]. Both are in agreement with the 
expected value of 6.0 x (1.000 ± 0.006) x 10^° cm“^s“^. 

The survival probability for electron neutrino from pp-reaction is 
Pee{Borex) = 0.64±0.12. This is the fourth energy range explored by Borex¬ 
ino, all the Borexino results on the electron neutrino survival probability are 
presented graphically in FigJS) 

Taking into account that Borexino and other experiments measurements 
are independent, the results can be combined: 

(ppp = (6.37 ±0.46) X 

The electron neutrino survival probability measured in all solar but 
Borexino experiment is Pee{other) = 0.56 ±0.08, combining it with Borexino 
one we obtain: 


11 









Reaction 

GS98 [31] 

AGS09 [35] 

Units 

cm ^ 

Measurement 

MeV/lu 

L 

xlO^e W-s-i 

pp 

5.98 ± 0.04 

6.03 ±0.04 

xlQio 

6.0 ±0.8 [11] 

6.6 ±0.7 [26] 
6.37 ±0.46 

13.10 

3.76 ± 0.28 

pep 

1.44 ±0.012 

1.47 ±0.012 

xl0“ 

1.6±0.3 1331 

11.92 

0.009 ± 0.002 

’Te 

5.00 ±0.07 

4.56 ±0.07 

xTo^ 

4.87 ± 0.24 E] 

12.60 

0.276 ± 0.014 

»B 

5.58 ±0.14 

4.59 ±0.14 

xl0“ 

5.25 ± 0.16 1361 

6.63 

1.57 ±0.05 

xlO-* 

hep 

8.0 ±2.4 

8.3 ±2.5 

>a(F 

< 231371 




2.96 ±0.14 

2.17±0.14 

Xio* 

GNO: 

< 7.4133] 



R>o 

2.23 ±0.15 

1.56 ±0.15 


i/p 

5.52 ±0.17 

3.40 ±0.16 



Table 2: The Standard Solar Model predictions (for high metallicity and low 
metallicity abnndances) and cnrrent experimental statns of the Solar nentrino 
flnxes measnrement. The limits are given for 90% C.L.. The corresponding 
energy release is calcnlated in the last colnmn 


Pee = 0.60 ± 0.07, 

well compatible with theoretical prediction of the MSW/LMA model 

U.OUl_Q Q42. 

All available measnrements of the solar nentrino flnxes are shown in 
Tabm The total energy prodnction in the solar reactions observed till now 
(by detecting corresponding nentrino flnxes) is 4.04 ± 0.28 W-s“^ in a good 
agreement with a total measnred Lq = 3.846 x 10^® W-s“^. There is not 
mnch space left for the nnknown energy sonrces, the 90% C.L. lower limit for 
the total energy prodnction (conservatively assnming zero contribntion from 
the not-observed reactions) is L^ot = 3.68 x 10^® W-s“h If one assnmes that 
snch an nnknown sonrce exists, its total power with 90% probability can’t 
exceed 0.15 x 10^® W-s“^. In other words no more than 4% of the total energy 
prodnction in the Snn is left for the nnknown energy sonrces, conhrming that 
the Snn shines dne to the thermonnclear fnsion reactions. 
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